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ABSTRACT 

This paper describes methods for reducing vibration 
in helicopter rotor blades by determining the opti- 
mum sizes and locations of tuning masses through 
formal mathematical optimization techniques. An 
optimization procedure is developed which employs 
the tuning masses and corresponding locations as 
design variables which are systematically changed 
to achieve low values of shear without a large mass 
penalty. The finite-element structural analysis of 
the blade and the optimization formulation require 
the development of discretized expressions for two 
performance parameters: the modal shaping parame- 

ter and the modal 9hear amplitude. Matrix expres- 
sions for both quantities and their sensitivity 
derivatives are derived in this paper. Three opti- 
mization strategies are developed and tested. The 
first is based on minimizing the modal shaping pa- 
rameter which indirectly reduces the modal shear 
amplitudes corresponding to each harmonic of air- 
load. The second strategy reduces these amplitudes 
directly and the third strategy reduces the shear 
as a function of time during a revolution of the 
blade. The first strategy works well for reducing 
the shear for one mode responding to a single har- 
monic of the airload but has been found in some 
cases to be ineffective for more than one mode. 

The second and third strategies give similar re- 
sults and show excellent reduction of the shear 
with a low mass penalty. 


LIST OF SYMBOLS 

A k Amplitude of the generalized force 

d istr ibut ion 

DAF*^ Dynamic amplification factor 

F- k Generalized force, ibf 

j } Force vector, Ibf 

f, f Q Objective function 

g, g 0 Constraint function 

I Set of modes included in procedure 

K Set of harmonics of airload included in 

procedure 

(Kj Elastic stiffness matrix, Ibf/in 

m Mass per unit length, Ibm/ in 

[M] Mass matrix, Ibm 

M Generalized mass, Ibm 
M n Tuning mass, Ibm 

MSP Modal shaping parameter 


Mtot 

Total tuning mass, Ibm 

NCON 

Number of constraints 

NDV 

Number of design variables 

NHARM 

Number of harmonics of airload included in 
procedure 

NMASS 

Number of tuning masses included in 
procedure 

NMODE 

Number of modes included in procedure 

NT I ME 

Number of critical time points in one 
revolution of blade 

<Uk 

Response of ith mode subject to k/rev 
harmonic of airload 

KJ 

Amplitude of the response of ith mode to 
k/rev harmonic of airload 

S 

Amplitude of the blade root vertical shear 
force, Ibf 

s ik 

Amplitude of the modal shear force 
corresponding to the ith mode and k/rev 
harmonic of airload, Ibf 

Sk 

Amplitude of the 3hear force associated 
with the k/rev harmonic of airload, Ibf 

9 

Time variation of blade root vertical 
shear force, Ibf 

8 ik 

Time variation of modal shear force 
corresponding to ith mode and k/rev 
harmonic of airload, Ibf 

c 

Time, seconds 

c m 

Critical time, seconds 

v i 

jth design variable 

3 max 

Maximum peak shear in one blade revolu- 
tion, Ibf 

X n 

Locations of the tuning masses along blade 
span, inches 

B,B k ,e. k 

Additional design variables 

h 

Damping coefficient 


Phase angle of k/rev harmonic 

{♦i> 

Eigenvector for ith mode 

Y ik 

Phase angle 

1 

Eigenvalue equal to the square of the ith 
natural frequency, rad/s 
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n Angular velocity of blade, rpm or rad/» The method is to formulate and solve an optimiza- 

tion problem in which the tuning masses and their 
Lower bound on frequency, rad/s locations are design variables that minimize a com- 

bination of blade shear and the added mass with 

W ui Upper bound on frequency, rad/s constraints on frequencies to avoid resonance. 


p Adjustable tolerance factor (see Eq . (37)) 

Av . Change in jth design variable 

Subscripts 

i Denotes ith mode 

j Denotes jth design variable 

k Denotes kth harmonic of airload 

Superscript 

T Denotes transpose of a matrix 


INTRODUCTION 

The current trend in engineering design of aircraft 
and spacecraft is to incorporate critical require- 
ments from alt pertinent disciplines analytically 
in an early phase of the design process to avoid 
the costly modification of a prototype after a 
problem has been detected (Ref. 1). Incorporation 
of vibration requirements in rotorcraft design is 
one example of this. In helicopter rotor blade and 
fuselage design, the need to increase ride comfort, 
stability, and fatigue life of structural compo- 
nents leads to stringent design constraints on 
vibration levels (Refs. 2-4). 

Vibration is transmitted from the blade to the 
fuselage primarily through a time-dependent shear 
force at the hub. Historically, frequency place- 
ment has been the principal technique for reducing 
rotor blade vibration (Refs. 5 and 6). Recently, 
it has been shown that minimization of the hub 
shear can be achieved through minimization of a 
modal shaping parameter (MSP) (Refs. 2 and 3). An 
associated technique, sometimes referred to as 
’'modal shaping" or "modal tailoring", alters the 
vibration mode shapes of the blades through mass 
and/or stiffness modificstion to make them less 
responsive to the airloads (Refs. 2, 7-9). A 
number of passive control techniques show promise 
for overall reduction of structural vibration in 
rotor blades. For example, pendulum absorbers 
(Ref. 10), active isolation devices (Ref. 11), 
additional damping (Refs. 4, 11-12), and vibra- 
tion absorbers which create anti-resonances 
(Refs. 13-14) have been demonstrated. Particularly 
effective is the strategic placement of tuning 
masses along the blade span to tailor mode shapes 
(Refs. 2-3, 5). What has been lacking in modal 
shaping and frequency placement methods is a 
systematic approach for predicting the best loca- 
tions for the tuning masses along the blade span. 

The purpose of this paper is to develop and demon- 
strate a method for optimally locating, as well as 
sizing, tuning masses to reduce vibration using 
formal mathematical optimization techniques. The 
design goal is to find the best combination of 
tuning masses and their locations to minimize blade 
root vertical shear without a large mass penalty. 


The optimization procedure includes a finite- 
element vibration analysis (Ref. 15) of a rotor 
blade in combination with a general-purpose optimi- 
zation code (Ref. 16). An explicit, approximate 
analysis (Ref. 17) of the blade vibration behavior 
is used to avoid the high computational cost of 
repeating the finite-element analysis for every 
blade modification. Sensitivity derivatives of the 
structural behavior with respect to the design 
variables are required in the optimization proce- 
dure. The research described herein includes de- 
velopment of the finite-element formulation of the 
analytical sensitivity derivatives of Che MSP and 
the blade root vertical shear. 

Three alternate optimization strategies are devel- 
oped and demonstrated in the paper. The first is 
based on minimizing the modal shaping parameter, 
thus reducing amplitudes of the modal shear for a 
single mode and single harmonic of the airloading. 
The second strategy directly reduces the shear am- 
plitudes corresponding to several harmonics for 
several modes, and the third strategy reduces the 
total shear as a function of time during a revolu- 
tion of the blade. Results are shown in which the 
above strategies are applied to a rotor blade de- 
sign considering single mode/single harmonic air- 
load, as well as multiple mode/multiple harmonic 
airload cases. 

ROTOR BLADE DYNAMIC ANALYSIS CONSIDERATIONS 

Calculation of the natural vibration mode shapes 
and frequencies, and the steady-state harmonic 
response are the fundamental analysis steps in the 
optimization procedure to be described. The pur- 
pose of this section of the paper is to outline the 
analytical basis and modeling conventions for the 
calculations . 

The rotor blade is modeled as a pinned-free beam 
undergoing lateral vibration normal to the plane of 
the rotor disk. This "flapwise" motion is in gen- 
eral accompanied by in-plane (edgewise) motion as 
well as torsion, but these are not included in this 
work. The beam is assumed to rotate at a constant 
speed (the rotor speed) about an axis which passes 
through the pinned end of the beam and is normal to 
the rotor disk. The effects of blade rotation are 
included through the centrifugal stiffness terms 
and differential stiffness terms in the equation of 
motion (see appendix). An additional term due to 
the Coriolis acceleration generally occurring in 
the equations of motion of rotating structures is 
not necessary in the model used herein. Also, 
damping is neglected in the calculation of frequen- 
cies and mode shapes and although it is often in- 
cluded in response calculations as modal damping, 
damping is neglected in all calculations herein. 
Based on the preceding, the governing equation for 
the finite-element modeled rotor blade is 

MX ♦ KX - F 
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where M if the mas* matrix 

K is the stiffness matrix 

F in the applied force vector 

Solution of the governing equation is described in 
more detail in the appendix, along with a discus- 
sion of how the centrifugal effects are included. 

The time-dependent forces acting on the rotor blade 
are associated with the lift and drag forces gener- 
ated by the air flow passing through the rotor 
dish. These loads are spatially distributed along 
the blade and vary with time sinusoidally at 
frequencies which are integer multiples (harmonics) 
of the rotor speed. It is customary in rotorcraft 
dynamics to use the notation N/ rev or N per rev 
when referring to frequencies or loading at N times 
the rotor speed. In rotor blade dynamic response, 
it is generally assumed that the loads are speci- 
fied in terms of distributions and phases for all 
appropriate harmonics. For an N-bladed rotor, the 
most critical frequencies of load and response are 
Nft and ( Nl t l)ft. Thus the most important harmonics 
are N and N t 1 . 

In the work described in this paper, the blade 
flapwise response is calculated by modal super- 
position. An outline of the derivation is given in 
the appendix, mainly for the purpose of establish- 
ing the terminology, but also for completeness. 

The flapwise modes which have the largest contribu- 
tion to vertical shear transmitted from the blade 
to the fuselage are the first and second elastic 
modes. Therefore these two modes are included in 
the example problems described in subsequent sec- 
t ions of the paper. 


DERIVATION OF PERFORMANCE PARAMETERS AND 
THEIR SENSITIVITY DERIVATIVES 

The degree to which the design is optimized is 
measured by two performance parameters: the modal 

shaping parameter (MSP) and the amplitude of the 
blade root vertical shear S. This section of the 
paper contains derivations of these performance 
parameters and their sensitivity derivatives. The 
roles they play in the optimization will be dis- 
cussed in a subsequent section of the paper. 


where is the steady state modal response 

(derived in appendix) 

u* is the ith natural frequency 

^ is the ith eigenvector 

and m is the mass per unit length 

For a discrete ( finite-element modeled) structure, 
the analogous expression is 

s ik ' q ik “i<«> T i*nv (2) 

where (M] is the diagonal mass matrix and ( U > is a 
Selection vector that extracts the appropriate com- 
ponents from the eigenvector. For example, for 
modal shaping of the flapwise modes, {U} contains 
I .0 in the rows correspond ing to the flap degree of 
freedom and 0.0 elsewhere. Using Equation (A8) 
from the appendix for the time dependence of 
and letting 

a. • <^{U} T [M1{4>.} d) 

gives 

s ik ■ s ik sinUne * V (4) 

where S ik • (5) 

ft is the angular velocity of the rotor 
blade 

t is time 

and represents a phase difference between the 
components of the response. This phase angle ac- 
counts for the lag between the modes, as well as 
the harmonics of the force. For cases in which 
there is no damping, y* k reduces to which is the 

phase angle for the kth harmonic (see appendix for 
more details of y- k ) . 

The time variation of the total shear results from 
summing over the modes and harmonics 


Derivation of MSP and Shear 

Reference 2 derives an expression for the vertical 
shear contribution, S£ k , from the ith flapwise mode 
due to the kth harmonic of the airload for a dis- 
tributed parameter model of a rotor blade. This 
paper develops its finite-element counterpart. 
Readers who are not familiar with Reference 2 may 
find the appendix helpful. 

From Equation (II) of Reference 2, the expression 
for s ](c is 


a . . * q . . • ok / . d x 

tk M ik t J i 


(l) 


s(t) « l l S. k s i n ( k !5t * Y ik ) (6) 

k i 


Substituting Equation (A9) from the appendix for 
[ q . k | and Equations (3) and (5) 


ik 




M. 

t 


• DAF.. • F.. 
ik ik 


( 7 ) 


where for the case of no damping (5*0) the 
dynamic amplification factor DAF^ is 
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DAF.. - 
ik 


(«./kOr 


[(ui./kfl) 2 - l] 


( 8 ) 


M. is Che generalized aata 

F ik * 8 the force, (Equation (A10)) 

written as 


F u 


V F /<V 


(9) 


A^ is the amplitude of the force and {F^} is the 
force distribution for the kth harmonic. Using 
Equations (7)-(9) and rearranging terms leads to 
the expression 


S 


ik 


(U) T [M|U.) . (F k } T {9.) 
M. 


DAF ik • \ 


(10) 


The total contribution from all the modes to the 
kth harmonic of shear is 

s k ■ l s u (l 


Following Taylor (Ref. 2), the first term of Equa- 
tion (10) is identified as the modal shaping 
parameter (MSP). Thus 

{U} T [M|{,.} . {F lt > T {*. } 

MS, 1 5 l - (12) 

lk M. 

I 

Equation (12) shows that an MSP exists for each 
mode shape and load case (harmonic). Since the MSP 
is a function of the mode shape, it is possible to 
reduce its value by tailoring the mode shape. For 
example, if the mode shape is made orthogonal to 
the force distribution, then the value of the MSP 
is zero. As seen from Equation (10), can be 

reduced by reducing MSP ik while limiting the size 
of DAF* The size of DAF ik is controlled by 
placing constraints on the natural frequencies as 
suggested in Reference 5. Equations (6), (10), 
(11), and (12) are the basis for the optimization 
strategies which will be discussed in a later 
section of the paper. 


Derivation of Sensitivity Derivatives of 
Performance Parameters 

Derivative of MSPifc - Beginning with Equation (12) 
the derivative of the MSP can be calculated as: 


r [S-]“i>) • <V T <*1 

■$} 


♦ {U}‘ 


♦ {UriMH*.} • (F 


j - l , 2 , • • • , NDV M3) 


where NDV is the number of design variables. In 
Equation (13) {3${/3vj} is the derivative of the 
ith eigenvector with respect to the jth design 
variable which is calculated analytically by 
Nelson’s method (Ref. 18) and (3M/3vj) is the 
derivative of the mass matrix which is calculated 
by finite differences. For the design variables 
used (masses and locations) the finite difference 
derivatives are exact. 

Derivative of - Using Equations (10) and (12) 
we can write the amplitude of the shear as 


S ik - MSP ik . DAF., 


(14) 


and then the derivative is 


3S .. 3MSP.. 

ik ik _ A 

-r * — r • DAF., • A. 

3v . 3v . ik k 

J J 


3 DAF . 


♦ MSP., 
ik 


ik 




• A, 


3MSP. 


where 


ik . 


(15) 


3v . 
J 


is taken from Equation (13) and the 


derivative of the dynamic amplification factor 
(DAF) of Equation (8) is 


3DAF., 

ik 

3v . 

J 


a(<# 2 .) 


(kor[(w./kn ) 2 - i] z 3v j 


(16) 


2 

In Equation (16) w. is the ith eigenvalue. The 
1 

derivative can be calculated analytically with 

the following equation (Ref. 19): 


. {, jTfaiKi . 2 IIM]] ( , J 

3v. l V 3v . i 3v . 

j L j j J 


(17) 


where the derivative of the stiffness matrix 

with respect to the jth design variable, is cal- 
culated by finite differences. 

Equations (13) and (15) are the expressions for the 
derivatives of the performance parameters which are 
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used in the sensitivity calculations for the three 
different optimization strategies. 


OPTIMIZATION FORMULATION 


Design Coa l 

The design goal ia to find the optimum combina- 
tion of tuning masse* ^ and their locations X n 
(fig. 1) to minimize blade root vertical shear 
while avoiding an excessive mass penalty. The 
method is to formulate and solve an optimization 
problem in which the tuning masses and locations 
are design variables that minimize the objective 
function which is a combination of a measure of 
vertical shear and the added mass. Additionally, 
constraints are placed on the frequencies to avoid 
resonance. It is noted that because this optimi- 
zation formulation involves mass as the objective 
function and frequencies and harmonic response as 
constraints, it falls into the general category of 
optimization problems discussed in Reference 20. 

In that reference, the complication of a disjoint 
design space was identified as occurring - the dis- 
jointedness being associated with noncontiguous 
regions of the design space on either side of res- 
onance points. It turns out that in the current 
formulation, the disjoint design space problem was 
avoided because in rotor blades, frequency changes 
due to varying tuning masses tend to be 9mall. 

Also, even when the disjoint problem occurs it may 
be dealt with by generating several designs from 
different starting points - a technique which is 
generally used when the presence of multiple local 
minima is suspected. 

Three optimization strategies will be described in 
this section of the paper. In each method, addi- 
tional design variables (8) are U9ed. They facili- 
tate the trade off between desired performance and 
excessive mass penalty. A general-purpose con- 
strained optimization program, CONMIN (Ref. 16), is 
used. CONMIN requires derivatives of the objective 
function and constraints. All derivatives are ob- 
tained analytically as will be shown in this 
section. 


Strategy I 

The objective function is 


These constraints express the requirement that the 
MSP^ be less than 0 ik . The convention is that 
for g < 0 the constraint is satisfied and violated 
otherwise. Large values of 0£ k make the con- 
straints easy to satisfy, but cause a large objec- 
tive function. Conversely, small values of 0^ 
result in a small objective function, but make the 
constraints more difficult to satisfy. The opti- 
mizer, therefore, will tend toward designs with the 
lowest possible values of 8' k and therefore low 
values of MSP^. Additional constraints include 
upper and lower bounds on the frequencies to avoid 
resonance 

»? . < w? < to 2 . (21) 

l l t u i 

or 

g » J./J. - 1 < 0 

and (22) 

g » 1 - < 0 


The required derivatives of the objective function 
and constraints are obtained by differentiating 
Equations (18), (20), and (22). The derivative of 
the objective function is 


0 


*n 


II 

dv 


■ • if •<* 

NMASS 


(23) 


3 ik 


Using the dimensionless form of the constraint 
function (Eq. (20)), the derivatives of the 
constraints are expressed as 


is_ 

dv . 
J 


r 2 

-MSP.. /87. 

ik ik 


i 3MSP ik 


8 U 


3v. 


v j “ 8 U 

(24) 

v j " X n M„ 


l l 

k i 


ik , 


NMASS 

l 

n* 1 


where Equation (13) is used for calculating the 
(Ig) derivative of MSP^ k . The constraints on the fre- 

quencies also need to be differentiated. Taking 
the derivative of Equation (22) gives 


where i is an element of the set of included modes 
(I) and k is an element of the set of included har- 
monics of airload (K) , \ are the tuning masses, 
and e- k are the additional design variables which 
also appear in the constraints defined below: 


MSP.. 

ik 


< 8. 


ik 


i C I 
kf K 


( 19 ) 


or in standard dimensionless form: 

g . Msp ik /a. k - i < o 


( 20 ) 


ro 


|s- 

3v . 
J 


9 ^ 2 

l 3 “i . 1 a “i 

_ __ and - 


u> 


J 


li 


3v . 
J 


0 ik 

( 25 ) 


Equation (17) is used for calculating the 

eigenvalue derivative, # 

J 
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m * 1,2 


NTIME 


( 32) 


Strategy II 

In this strategy the constraints are placed on the 
harmonic amplitudes, (see Eq. (11)). The 
objective function is 

NMASS 

' -(■ *!•<*) J, M - 

The additional design variables 8^ play a role 
similar to 8^ in the first strategy. Here the 
constraints are written as 

S k < 0 k k i K (27) 

or 

g - V 3 k - 1 < 0 r23) 


This strategy also employs upper and lower bound 
constraints on the frequencies (Eqs. (21) and 
( 22 )). 

Similar to strategy l, the derivative of the objec- 
tive function is obtained by differentiating 
Equa t ion (26). 


. <: 

3v i 


1 ♦ I 8. 


NMASS 

l M n 
n»l 


v « X 
J " 

v j ■ M n 


v j ’ 8 k 


(29) 


Differentiating Equation (28) gives the derivatives 
of the constraints 


\S 3v i 


Vj - 8 k 

(30) 

v j * x n ° r \ 


k c ik ik . . 

where * l — , and — — is given in Equa- 

J i j j 

t ion (15). Equations (25) and (17) give the 
derivatives of the frequency constraints. 


Strategy III 

In this formulation, the shear as a function of 
time is minimized by constraining all the peak 
values which occur during a revolution of the 
blade. These are called critical point constraints 
(Ref. 21). The objective function is 

NMASS 

f » (l ♦ 8] l M (3l) 

n-1 


and the constraints are 

s(t ) < 0 
m 

or 

g - s( t ) / 8 - 1 < 0 (33) 

m 


which require that the values of the shear at each 
time t m be less than the value of 0; and again 0 is 
minimized because of its role in the objective 
function. As shown in Equation (6), s C t m ) is the 
shear at time t m where t m represents a time at 
which a peak occurs in s(t). 


The peak values are identified as follows: A value 

of NTIME is specified corresponding to the maximum 
number of peaks in the function s( t) during a revo- 
lution of the blade. The peak values are identi- 
fied by examining the shear as a function of time 
(Ref. 21). The constraints are placed on these 
NTIME values of shear to force the peaks to be as 
small as possible. This procedure does not require 
that the locations of the peaks be constant 
throughout the optimization process since the 
search for the peaks occurs each time the analysis 
is performed for a new set of design variables. 
Additional constraints are again placed on the 
upper and lower values of the frequencies 
(Eqs. (21) and (22)). Again the objective function 
for strategy III yields a similar derivative as in 
the two previous cases 


3f 

377 


(nmass 


. 


(\ + & 


0 


v j - *n 
v j - *n 


(34) 


Differentiating Equation (33) with respect to the 
design variables gives* 


-s(t m )/8 


2 



3s( t ) 
m 

3v . 

J 


0 


X n ° r M n 


(35) 


where s( t m ) is obtained from Equation (6) and 
m * 1 , 2 , . . . , NTIME. The derivative of the shear 
3s(t ) 

m 

— r at time t m is 

dv . " 


3s(t ) 
m 


3v. 

J 


J * X. 

i l S in( ktlt 

k i j 


m 


V 


where the derivative 
Equat ion (15). 


3S.. 
ik 

3v . 

J 


is given by 


(36) 


♦The critical point t m is also a function of the 
design variables; however, as pointed out in 
Reference 21 this has no effect on the derivative 
in Equation (35). 
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KS-Func t ion 


It cm be seen from the optimization formulations 
that as we begin to consider the contributions of 
larger numbers of modes and load cases to the 
shear, or as the number of critical points in 
strategy III increases, the number of constraints 
will becosie large. This can cause slower conver- 
gence of Che optimization process. An envelope 
constraint function, denoted the Kre isselmeier- 
Steinhauser (KS) function (Ref. 22), is used to 
substitute a single constraint function for a large 
number of constraints. The KS-function is defined 
as 


loop consists of the optimization program, CONMIN 
(Ref. 16) and an approximate analysis for calculat- 
ing the objective function and the constraints (see 
Ref. 17). The approximate equations are 


NDV 


e - e ♦ l 

n *• 


3f 


Av . 

. , jv . j 

NDV 

g * g + l Av . 

6 6 o . L . 3v . j 

J*l J 


OR) 


(40) 


KS - ~ in 
P 


/ NCON 

l 

\i-l 



(37) 


where the g^ are the actual constraints, NCON is 
the number of constraints, and p controls the dis- 
tance between the KS-function and the actual con- 
straint boundary (Fig. 2). For small values of p, 
the KS-function is very smooth and also a very 
conservative estimate of the constraint viola- 
tions. As p increases, the KS-function moves 
closer to the discontinuous function, max (g:). 

i 1 

Typically, the initial value of p is small and is 
increased as convergence is approached. 

The optimizer requires the derivative of the KS- 
function. Taking the derivative of Equation (37) 
g ives 


3KS 

3v . 
J 


NCON 

l 

i»l 


p *i 3 *i 


3V j 


NCON pg. 

I e 1 

i-l 


(38) 


where 



are derivatives of the constraints. 


OPTIMIZATION PROCEDURE 


The sequence of operations in the optimization 
procedure is illustrated in Figure 3. The overall 
procedure consists of two nested loops. Each pass 
through the outer loop is referred to as a cycle 
which involves a full analysis and a sensitivity 
calculation. The first step is to generate the 
f i n L te-e lement structural model of the beam, ex- 
cluding the values of tuning masses. The design 
variables (masses and locations) are used to allo- 
cate the masses to the appropriate grid points of 
Che model. Specifically, the masses are divided 
between the two grid points adjacent to each X n by 
prorating according to the distance from each. 

Next, the masses are inserted into the model, the 
vibration analysis is performed, and the MSP^ and 
shear amplitudes are calculated for NMODE 

number of modes responding to NHARM number of 
harmonics. The sensitivity analysis includes 
calculating the vibration mode shape derivatives by 
Nelson's method and then calculating the deriva- 
tives cf the objective function and constraints 
using Equations (23)— (25) for strategy I, Equations 
(25) and (29)-(30) for strategy II, and Equations 
(25) and (34)-(36) for strategy III. The inner 


These equations give Che change in the objective 
function from f Q to f and the change in a con- 
straint from g 0 to g corresponding to a change In 
design variable Av j . To assure that the linear 
approximation in Equations (39) and (40) are valid, 
the size of Avj is limited to 10 percent of vj ♦ 

Use of these approximations saves computational 
time and effort in the inner loop where many eval- 
uations of the objective function and constraints 
are required. Development of these and other tech- 
niques and demonstration of their benefits are 
described in Reference 23. It was observed that 
the KS-function of the linear approximation to g 
wa9 a very accurate approximation to the KS- 
function of the exact g. Once the inner loop 
iterations have converged the next cycle of the 
outer loop begins, using the current design vari- 
ables as the new values of the lumped masses and 
their locations. These masses are Chen inserted 
into the structural model and the process continues 
until convergence of the outer loop is achieved. 


EXAMPLE PROBLEM 


The example problem is a beam representation of an 
articulated rotor blade developed in Reference 5 
and shown in Figure 4, The beam is 193 inches long 
with a hinged end condition and is modeled by 10 
finite elements of equal length. The model con- 
tains both structural mass and lumped (non- 
structural) masses. The beam has a box cross 
section as shown in Figure 4b and the material 
properties and cross-sectional dimensions are sum- 
marized in table 1. Three lumped masses are to be 
placed along the length of the beam. The values of 
the masses and their locations are the design vari- 
ables and their initial values, shown in table l, 
are from the blade in Reference 5. 

Figure 5 shows the distributions (from Ref. 2) and 
phase angles (from Ref. 24) of the airloads used. 
They are input as tabulated values of distributed 
forces (i.e., force per unit length) into the 
finite-element analysis (Ref. 15). The forces 
needed in Equation (9) are calculated as 
consistent nodal forces for the finite-element 
model. These airloads represent the 3, 4, and 5 
per rev lifting airloads that are typical of a 
four-bladed articulated rotor system. These 
harmonics were chosen since they are the prime 
contributors to the vibration of a four-bladed 
rotor system. In this work, the first and second 
elastic flapping modes are included since they are 
prime contributors to the vertical vibration 
transmitted by the rotor blades to the fuselage. 
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RESULTS AND DISCUSSION 

The following sections discuss results obtained for 
each of the three optimisation strategies applied 
to the example problem. The test cases include ( 1 ) 
a single mode responding to a single harmonic of 
the airload; ( 2 ) two modes responding to a single 
harmonic of airload; and (3) two modes responding 
to three harmonics of the airload. 

Results for Strategy I 

Figures 6 a and 6 b show the initial and final de- 
signs for the first test case using strategy l for 
minimizing Sj^, the shear for the first elastic 
flapping mode and the 4/rev airload by proper 
placement of three tuning masses. The 4/rev air- 
load is concentrated at the tip of the blade and in 
order to shape the mode to be insensitive to the 
airload all three tuning masses were moved to the 
tip of the blade. Figure 6 c shows a sketch of the 
mode shape before and after optimization. It is 
this change in the mode shape that reduces the 
value of the MSP^ 99 percent and the corresponding 
shear by a similar amount. Table 2 summarizes 

the initial and final designs. It is noted that 
the large change in MSP and shear is accompanied by 
an almost zero change in frequency. 

This method proved to be very useful when working 
on the response of one mode corresponding to one 
airload; however, when this method was applied to 
more than one mode it was not always effective. 

The method did reduce the values of the MSP's as 
required, but low values of the MSP's did not 
necessarily give low values for the shear, Sy The 
reason for this is that the various MSP's (and, 
therefore, the corresponding contributions to the 
shear, S i(c ) may have different signs. When these 
shear contributions are added together (see Eq . 
(ID), small values for the individual contri- 
butions do not always minimize the shear Sy 
unless the necessary cancelling of equal and 
opposite terms occurs. 


Results for Strategy II 

In strategy II, the method is to reduce the values 
of the shear contributions, S^ summed over the 
modes (see Eq . (ID). This eliminates the problem 
in strategy I since constraints placed on the sums 
of the modal shears encourage the desirable can- 
cellation effects. The first step in demonstrating 
this method was to validate it for the 1 mode /1 
load case. The design for strategy II is essen- 
tially identical to that of Figure 6 and table 2 
gives the initial and final designs of the blade. 
The next step was to apply this strategy to a 2 
mode/l load case. The response of the first and 
second elastic flapping modes corresponding to the 
4 /rev airload is minimized. Table 3 summarizes the 
initial and final designs. T^e initial shear 
including both modes is -34.68 lb f which is reduced 
by the optimization process to -.01 lb f with an 
accompanying decrease in the added mass. 

Inspection of the initial and final values of the 
MSP's (table 3) shows that the magnitudes of the 
MSP's became larger. This helps to explain the 
Lack of success of strategy I for this problem. In 
fact, the modal shears Sj^ and S 2 ^ are equal and 


opposite thus combining together to produce a near 
zero value of total shear. Figures 7 a and 7b show 
the initial and final masses and their locations 
and Figure 7c shows the change in the shapes of the 
first and second elastic flapwise modes. 

Strategy II was next applied to a case of two modes 
responding to the 3, 4, and 5 per rev harmonics of 
airioading (see Fig. 5). Table 4 shows the initial 
and final results where the amplitudes of the 
shears, due to the two flapwise modes, have been 
reduced significantly with only a 9 Ibm increase in 
tuning mass. For example, S 5 , the shear force 
associated with the 5/rev harmonic, wa 9 reduced 
from -39.48 Ibf to -.162 tbf. Figure 8 gives a 
time history of the shear during one revolution of 
the blade before and after optimization. 

It is clear from the figure that reducing the 
amplitudes of the harmonic shears results in a 
large reduction of the total shear throughout i 
revolution of the blade. 


Results for Strategy III 

The third strategy was applied to the previous test 
case of two modes responding to three harmonics. 
Figure 9 shows graphs of the shear s( t) plotted as 
a function of time and azimuth for a revolution of 
the blade for the initial and final designs. The 
peaks on the initial curve have been reduced 
dramatically. For example, the maximum peak s max 
for the initial design is -78.00 Ibf and for the 
final design the maximum peak is -.576 Ibf. The 
extreme right column in table 4 gives details of 
the final designs from strategy III and indicates a 
large payoff for a relatively small Increase in 
added mass. 


Comparison of Strategies II and III 

As seen from table 4, strategy III gave a design 
comparable to strategy II. In both cases, there 
was a significant reduction in the total shear. 
Results from strategies II and III are compared in 
table 5 , in terms of the peak values of s( t) at the 
critical points. The final results were very close 
with strategy III producing slightly lower values 
for some of the peaks and strategy II producing 
lower values of others. Overall, strategy lit was 
slightly better in minimizing the peak shear. 

Figure 10 shows the shears plotted as functions of 
time from each strategy. Strategy III was a some- 
what more complicated approach than strategy II and 
was more cumbersome to implement. However, the 
greater degree of rigor in strategy III makes the 
slightly greater effort worth the investment. 


CONCLUDING REMARKS 

This paper described methods for systematically 
locating, as well as sizing, tuning masses to 
reduce vibration in helicopter rotor blades using 
formal mathematical optimization techniques. The 
problem was to find the optimum combination of 
tuning masses and their locations to reduce ver- 
tical shear without a large mass penalty. The 
methods embodied optimization procedures in which 
tuning masses and their locations were design 
variables whose values minimize a combination of 
shear and added mass. The finite-element 
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structural analysis of the blade and the optimiza- 
tion formulations required the development of dis- 
cretized expressions for two performance parame- 
ters: the modal shaping parameter, and the ampli- 

tude of the blade root vertical shear. Matrix 
expressions for both quantities and their sensitiv- 
ity derivatives were derived in this paper. The 
mechanism in the optimization for reducing the ver- 
tical blade shear was through "modal shaping" by 
placing the tuning masses at strategic locations 
along the blade. 

Three optimization strategies were developed: the 

first was based on minimizing the modal shaping 
parameter, thus indirectly reducing the amplitudes 
of the modal shear for each harmonic; the second 
reduced the shear amplitudes directly; and the 
third reduced the total shear as a function of time 
during a revolution of the blade. 

Strategy I worked well for reducing the shear for 
one mode responding to one harmonic of the load, 
but was ineffective for multiple modes. This was 
due to inability of the method to take advantage of 
sign differences between the contributions for 
different modes and harmonics. Strategy II worked 
extremely well for the 1 mode/1 load case, as well 
as multiple mode/multiple loads. Strategy III gave 
excellent results; that is, the peak shear was re- 
duced significantly without a large mass penalty. 
Strategies II and III gave esentially the same 
results for a 2 mode/3 load case. Strategy II is 
slightly easier to implement but the fact that 
strategy III is a more rigorous approach make9 it a 
preferred choice overall. 


APPENDIX 

D ERIVATION OF GENERAL FORM OF STEADY-STATE 
RESPONSE OF A HELICOPTER ROTOR BLADE 


where K £ is the linear elastic stiffness matrix 

is a centrifugal stiffness matrix which 
contains products of masses and angular 
veloc ity 

K d is the differential stiffness matrix and 
contains stresses associated with exten- 
sion of the beam due to centrifugal 
forces 

Details of the derivations and explicit form9 for 
K c and K d may be found in References 25 and 26. 

Express the response vector X as a modal expansion 
such that 

NMODE 

X * F $ • q • (A3) 

l* l 

where is the ith vibration mode of the 

structure 

q £ is the generalized response vector 
The modes are normalized such that 

*7 M*. * M. (A4) 

T i T i i 

Equation (A4) is the usual definition of as the 
generalized mass associated with the ith mode. 
Combining Equations (Al), (A3), (A4) gives 

q. ♦ 2£.ui.q. + w?q . * — F. (A5) 


The purpose of this appendix is to develop the 
general form of the steady-state modal response q^ 
for the ith mode and kth harmonic of loading. 

This derivation is based on that of Reference 2 and 
is included here to establish the notational con- 
ventions used in the main body of the paper in the 
expressions for the modal shaping parameter (MSP) 
and blade root vertical shear (S). 

The governing matrix equation for vibration re- 
sponse; of a finite-element modeled structure is 


where is the natural frequency of the ith mode 

F^ is the generalized force for ith mode 

is modal damping coefficient (see 
' Ref. 27) 

Equation (A5) is identical to Equation (1) of 
Taylor (Ref. 2). The generalized force vector F ^ 
is a superposition of contributions from all 
harmonics . 


MX ♦ CX * KX - F 
where M is the mass matrix 


(Al) 


l F ik cos k^t * l F^ ks sin kfit (A6) 
k*0 1 C k* 1 


C is the damping matrix 

K i9 the stiffness matrix 

X is the vector of displacements and 
rotations 

F is the applied force vector 

The stiffness matrix K for a rotor blade has the 
form 


K * k E * K C * k D 


( A2 ) 


From Equation (A5) and (A6) it follows that the 
steady state modal response amplitudes q^ are of 
the form 


q i ” I q ikc cos kJlt * I q ik* sin kQt (A7) 
k-0 k-1 


■ l lq.^1 ♦ V (A8) 

k-0 
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The quantity | q - . | is the amplitude of the re- 
sponse. The defining equation for | q ^ ^ | is 




•u' 




T72 


(A9) 


where 


,2 t 1/2 


F., - [F Z . • FT. ] 
Lk 1 iks ikc J 


Y ik ’ hk + V 


( A10) 


where 


-1 


hk * tan 


25 i w i /kn 


(!i) . 

,\ kfl/ 


(All) 


Equation (AiO) along with Equation (All) show that 
the phase angle y^ is composed of the sum of two 
contributions. The first is the usual phase 

due to damping of the ith mode and the second is 
due to the phase lag in the kth harmonic of load. 
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Table 1. Details of finite-element model 
of rotor blade. 



(a) Material properties and 
cross-sectional dimensions 
( see F igure 4) 


Element 

No. 

E p b h 

(psi) (lb/in^) (in) (in) 

t d 

( in) ( in) 

1 

2-10 

0.490 * 10 7 .07 3.75 2.50 

0.585 x 10 7 .07 3.75 2.50 

o o 
o o 

CO 00 

(b) 

Initial values of design variables 
(see Figure 1) 

( in) 

X 2 X 3 Mj M 2 

( in) ( in) ( Ibm) ( Ibm) 

M 3 

( Ibm) 

135.10 

154.40 173.70 5.21 8.55 

6.60 


(c) Blade characteristics 


Length 
( in) 

Structural Nonstruc tura l 

mass (Ibm) mass (Ibm) 

a 

( rpm) 

193 

87.82 122.93 

425 

Table 2. Initial and final designs using 
strategies I and II for l mode/i load problem. 
(First elastic flapwise mode at 4/rev) 


Initial Final 



Strategy I Strategy II 


X t ( in) 

135.10 

193.00 

193.00 

X 2 ( in) 

154.40 

193.00 

193.00 

X-j ( in) 

173.70 

193.00 

193.00 

Mj (Ibm) 

5.21 

6.30 

6.70 

M-> ( Ibm) 

6.55 

10.53 

10.33 

( Ibm) 

6.60 

10.60 

10.40 

Mtot ( Ibm) 

18.36 

27.43 

27.43 

MSP,, 

-11.75 

-0.008 

-0.010 

O) 
c r 

— h 

-58.50 

-0.040 

-0.051 

w, (per rev) 

2.69 

2.70 

2.70 


Table 3. Initial and final designs using 
strategy II for 2 mode/1 load problem. 

(First and second elastic flapwise modes at 4/rev) 


Initial 

Final 

(in) 

135.10 

135.02 

X 2 ( in) 

154.40 

146.38 

X 3 (in) 

173.70 

157.72 

Mj ( Ibm) 

5.21 

7.15 

M 2 (Ibm) 

6.55 

5.24 

Mj ( Ibm) 

6.60 

0.97 

Mtot ^ l^ m ) 

18.36 

13.36 

s 4 (Ibf) 

-34.68 

-0.01 

S, 4 (Ibf) 

-58.50 

-66.91 

S 24 (Ibf) 

23.82 

66.90 

msp 14 

-11.75 

-14.61 

MSP 24 

1.03 

2.24 

U] (per rev) 

2.69 

2.63 

id, ( per rev) 

4.65 

4.47 

Table 4. Initial 

and final 

designs using 


strategies II and lit for 2 mode/3 load problem. 
(First and second elastic flapwise modes it 
3, 4, and 5 per rev) 



Initial 

Final 


Strategy II 

Strategy III 

X t ( in) 

135.10 

150.95 

152.81 

x 2 (in) 

154.40 

154.44 

154.27 

X 3 (in) 

173.70 

154.40 

154.42 

Mj ( Ibm) 

5.21 

4.28 

7.93 

M 2 (Ibm) 

6.55 

12.40 

9.66 

M 3 (Ibm) 

6.60 

10.13 

9.26 

Mtot ( ibm) 

18.36 

26.81 

26.85 

S 3 (Ibf) 

-7.98 

-0.36 

- 

s 4 (ibf) 

-34.68 

0.30 

- 

s 5 (Ibf) 

-39.48 

-0.16 

- 

*max < lbf > 

-78.00 

- 

-0.58 

u»l (per rev) 

2.69 

2.81 

2.81 

u >2 ( per rev) 

4.65 

4.72 

4.72 


Table 5. Comparison of peak values of shear 
during one revolution of the blade using 
strategies It and III. 

Strategy It Strategy tit 


.<«„> ate) (ini t .(.) 


-.348 

.009 

.330 

.000 

.198 

.024 

-.186 

.014 

.162 

.044 

.234 

.042 

-.528 

.059 

-.390 

.059 

.756 

.075 

.510 

.075 

-.768 

.092 

-.576 

.090 

-.588 

.107 

.570 

.107 

-.402 

.121 

-.504 

.123 

.354 

.136 

.366 

.138 

.132 

.141 

- 

- 
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Figure l. Problem definition: selection of 

optimum locations for tuning masses 
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(b! Cross sectional detail of r otor blade 




Figure 3. Flow chart for optimization procedure 


Figure 4. Test problem 



Percent blade span 


(0 3 = 59°) 
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Figure 5. Amplitude distributions and phase 
angles of airloads 
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(a) Initial masses and locations 
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• 
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(b) Final masses and locations 



Figure 6. Initial and final designs for 1 mode/ 
l load case using strategy l 



t/T 

Figure 8. Time history of blade vertical root 

shear during one revolution. Initial 
and final designs for 2 mode/3 load 
case us ing strategy II 
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Initial 

Final 


135. 154. 174. uni Mode 1 
tai Initial masses and locations 


Initial shear. S 4 = -34.74 Ibf 


Final shear. $ 4 = ;0054 Ibf 



7 1 5.2 Q7 
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<h» Final masses and locations 
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* c> Initial and final mode shapes 


Figure 7. Initial and final designs for 2 mode/ 
l load case using strategy II 
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Figure 9. Time history of blade vertical root 

shear during one revolution. Initial 
and final designs for 2 mode/3 load 
case using strategy III 
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Figure 10. Comparison of final designs from 
strategies II and III for 2 mode/ 
3 load case 


14 



J|Q§£Li Report Documentation Page 

Soace VsrraMon 

1. Report No. 

NASA TM- 100562 
AVSCOM TM 88-B-003 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Optimal Placement of Tuning Masses 
for Vibration Reduction 

5. Report Oat. 

March 1988 

6. Performing Organization Code 

7. Author(s) 

Jocelyn I. Pritchard and 
Howard M. Adelman 

8. Performing Organization Report No. 

10. Work Unit No. 

505-63-51-10 

9. Performing Organization Name and Address 

NASA Langley Research Center, Hampton, VA 23665 and 
Aerost ructures Directorate, USAART A- AVSCOM, 

Hampton, VA 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration, 
Washington, DC 20546 and U.S. Army Aviation Systems 
Command, St. Louis, MO 63120-1798 

14. Sponsoring Agency Code 

15. Supplementary Notes 


vlocelyn I. Pritchard, Aeros tructures Directorate, USAARTA-AVSCOM, Langley Research 
Center, Hampton, VA. Howard M. Adelman, Langley Research Center, Hampton, VA. 
Presented at the 29th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and 
Materials Conference, April 18-20, 1988, Williamsburg, Virginia. 

16. Abstract Th e paper describes methods for reducing vibration in helicopter rotor 
blades by determining optimum sizes and locations of tuning masses through formal 
mathematical optimization techniques. An optimization procedure is developed which 
employs the tuning masses and corresponding locations as design variables which are 
systematically changed to achieve low values of shear without a large mass penalty. 
The finite-element structural analysis of the blade and the optimization formulation 
require development of discretized expressions for two performance parameters: 
modal shaping parameter and modal shear amplitude. Matrix expressions for both 
quantities and their sensitivity derivatives are derived. Three optimization 
strategies are developed and tested. The first is based on minimizing the modal 
shaping parameter which indirectly reduces the modal shear amplitudes corresponding 
to each harmonic of airload. The second strategy reduces these amplitudes directly, 
and the third strategy reduces the shear as a function of time during a revolution 
of the blade. The first strategy works well for reducing the shear for one mode 
responding to a single harmonic of the airload, but has been found in some cases to 
be ineffective for more than one mode. The second and third strategies give similar 


results and show excellent reduction of the shear with a low mass penalty. 


17. Key Words (Suggested by Author! s)) 

18. Distribution Statement 

Modal shaping 


Optimization 

Unclassified - Unlimited 

Rotor blade vibrations 

Subject Category 39 

Structural dynamics 


Mass Placement 


19. Security Classif. I of this report 1 
Unclassified 


| 20. Security Classif. (of this page) 
Unclassified 


21. No. of pages 

15 


NASA FORM 1626 OCT 86 

















